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ABSTRACT 

An  extension  of  the  theoretical  work  of  Lamb  and  Pekeris 
shows  that  a  prominent  phase  (Airy  Phase)  consisting  of  a  group  of 
waves  of  period  9“H  seconds,  travelling  across  the  ocean  with  a 
velocity  of  approximately  0*7  of  the  speed  of  sound  in  water  should 
he  present  on  the  seismograms  of  shallow  focus  submarine  earth¬ 
quakes.  This  arrival  corresponds  to  normal  mode  propagation  at  a 
stationary  value  of  group  velocity  through  the  acoustic  system 
consisting  of  the  ocean  and  the  rigid  ocean  bottom. 

A  phase  fitting  this  description  has  been  observed  on 
the  Milne-Shaw  instruments  at  Bermuda  for  a  series  of  Dominican 
Republic  shocks.  The  Wenner  seismographs  at  Huancayo  recorded  the 
Airy  phase  of  a  shock  southwest  of  the  Galapagos.  Islands. 

The  study  of  the  Airy  phase  of  a  sufficient  number  of  sea 
quakes  would  provide  information  concerning  the  nature  of  the  ocean 
bottom  since -the  -latter  affects  both  the  periods  and  velocities  of 
the  Airy  phase. 


2. 


Introduction 

A  number  of  seismograms  of  earthquakes  occuring  at  sea 

were  examined  and  it  was  found  that  for  a  series  of  these  a  new 

phase  could  he  identified  consisting  of  a  group  of  waves  of  period 

8-9  seconds,  travelling  across  the  ocean  with  a  velocity  of  0*7- 

0.8  of  the  speed  of  sound  in  water.  It  is  believed  that  this 

arrival  corresponds  to  normal  mode  propagation  at  a  stationary 

value  of  group  velocity  through  the  acoustic  system  consisting  of 

the  ocean  and  the  rigid  ocean  bottom. 

Pekeris^  predicted  and  later  identified  this  phase  on 

seismograms  of  explosive  sounds  in  shallow  water  obtained  by 

2 

Worzel  and  Ewing.  In  his  theory  of  normal  mode  propagation  in 
two  liquid  layers,  Pekeris  named  this  arrival  the  Airy  phase  and 
showed  that  its  period  depended  on  the  depth  of  the  water  and  the 
velocity  contrast  between  the  water  and  the  bottom.  That  periods 
associated  with  a  stationary  value  of  group  velocity  would  pre¬ 
dominate  at  long  distances  had  been  recognized  earlier  by  Jeffreys^ 
h 

and  Stoneley. 

A  theory  for  the  origin  of  micro  seisms  in  which  micro- 

seism  groups  are  considered  to  be  Airy  phases  was  advanced  by 

5 

Press  and  Ewing.  This  theory  explains  the  mechanism  of  transmission 
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of  seismic  energy  from  a  storm  at  sea  and  correctly  predicts  that 
periods  of  three  to  nine  sec  will  predominate,  the  period  being 
determined  by  the  depth  of  water  and  the  elastic  constants  of  the 
bottom. 

Whether  the  source  of  seismic  energy  is  a  disturbance 
at  the  ocean  surface  produced  by  a  hurricane,  an  explosive  charge 
detonated  within  the  ocean,  or  an  earthquake  beneath  the  bottom, 
the  normal  mode  theory  indicates  that  a  prominent  Airy  phase  should 
be  observed  at  a  distant  point.  Figure  1  ,  the  seismogram  from  a 
shot  of  5  lbs.  in  water  57  feet  deep,  at  a  distance  of  4.5  nautical 
miles,  illustrates  the  predominance  of  the  Airy  phase  for  shallow 
water  explosions. 

Although  the  treatment  of  the  bottom  as  a  liquid  was 
adequate  for  the  shallow  water  data  to  which  Pekeris  applied  it, 
when  one  deals  with  submarine  earthquakes  in  the  deep  ocean  it 
becomes  necessary  to  consider  the  bottom  as  an  elastic  solid.  In 
this  paper,  the  Pekeris  theory  is  extended  to  include  the  case  of 
an  impulsive  point  source  of  compressional  waves  within  the  bottom 
which  is  taken  as  a  solid. 

Propagation  of  Elastic  Waves  in  the  Ocean  for  an 

Impulsive  Point  Source  beneath  the  Bottom 

The  analagous  problem  of  the  generation  of  Rayleigh 

waves  at  the  free  surface  of  a  semi-infinite  solid  for  a  source 

c 

within  the  solid  was  studied  by  Lamb;  Subsequent  investigations 

6.  Horace  Lamb,  On  the  Propagation  of  Tremors  over  the  Surface  of 
an  Elastic  Solid,  Phil.  Trans.,  Roy.  Soc.  Lon.,  Series  A, 
Vol.  203,  PP  1-^2,  1904 
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by  Stoneley^  and  Sezawa®  dealt  with  the  effect  of  the  ocean  on 
the  transmission  of  seismic  waves  across  the  ocean  bottom.  Stoneley 
considered  plane  waves  and  calculated  several  values  of  phase  and 
group  velocity  for  long  periods.  Sezawa  obtained  an  approximate 
solution  for  the  propagation  of  cylindrical  waves  having  great 
length  compared  to  the  depth  of  the  water.  Lee^  calculated  the 
amplitudes  of  Rayleigh  waves  for  plane  waves  transmitted  through 
a  homogeneous  medium  covered  by  a  superficial  layer,  assuming  various 
values  for  the  elastic  constants  of  the  two  media.  Scholte^  gave 
a  theory  for  the  combined  effects  of  gravity  and  compressibility  in 
a  layer  of  water  in  contact  with  an  elastic  solid  bottom.  Although 
Scholte’s  emphasis  was  not  on  propagation,  he  did  prepare  an  ex¬ 
cellent  tabulation  of  phase  velocities  for  the  first  made.  Press 
and  Ewing^  presented  complete  curves  of  phase  and  group  velocity 
for  the  first  and  second  modes  of  motion  in  a  liquid  layer  super¬ 
posed  on  a  semi-infinite  solid  bottom  and  noted  the  importance  of 
stationary  values  of  group  velocity  for  elastic  waves  transmitted 
across  the  ocean. 

In  the  present  paper  we  make  particular  use  of  the  methods 
and  results  of  Lamb  and  Pekeris  to  derive  the  normal  mode  solutions 
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pp  669-683,  1943 

11.  Frank  Press  and  Maurice  Ewing,  loc.  cit . 
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from  which  the  wave  motion  at  large  distances  from  the  source  can 
be  obtained. 


Consider  the  propagation  of  elastic  waves  through  the 
system  consisting  of  a  liquid  layer  of  thickness  H,  density  p ^ 
and  compressional  wave  velocity  v^,  superposed  upon  an  infinitely 
thick,  solid  bottom  of  density  whose  compressional  distortional 
wave  velocities  are  respectively  o^and^.  X/  is  the  compressi¬ 
bility  of  the  liquid  and  "X^ ,  ^  are  Lamfe’s  constants  for  the 
bottom,  and  are  related  to  the  elastic  wave  velocities  as  follows; 


The  Cartesian  coordinate  system  is  chosen  with  the  x,  y  axes  in 
the  free  surface  and  the  z  axis  vertically  downward.  Assuming 
symmetry  about  the  z  axis  we  have  r  =  J x^  +  y^  and  denote  the 
displacements  perpendicular  and  parallel  to  that  axis  by  q  and  w 
respectively.  The  subscripts  1  and  2  refer  to  the  liquid  layer 
and  bottom  respectively. 

Simple  Harmonic  Compressional  Point  Source  within  Bottom 


We  wish  to  determine  the  component  displacements  q  and 
w  due  to  a  point  source  of  compressional  waves  at  r  =  0,  z  =  d  «f  H. 
In  particular  we  will  seek  the  normal  mode  solutions  which  pre¬ 
dominate  at  large  distances  from  the  source. 

Following  the  procedure  of  Lamb,  we  assume  simple  harmonic 
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motion 


in>T 

e 


and  introduce  the  functions  (pfcz)  and  ^(r,  z)  defined 


by  the  equations: 

$,  -  c,^/°,r 
r  dlf./d* 


ur. 


%  -  J'/k/j'r  +c f'Pz/drJz 

<4  =  ^afdz  +  ^%/dzx  +  rfg*  <4 


The  vertical  stress  pzz  and  the  tangential  stress  pZT  can  be  ex¬ 
pressed  in  terms  of  q,  w,  (/? ,  (/y  and  the  elastic  constants  follow: 


P«-  XvV  -va./f&ur/d z  - 

Pzr=  ju  (  s  ifdz  +  2  ur/d  r) - 


(l) 

(«) 


It  is  required  that  the  functions  (0  and  IjJ  satisfy  the 


equations : 

(  Va+  If,  =  o  - - - (F) 

( va  +  uy°<i)iPt  =  o  . . -  -  (a) 

( va  +  u)a/4aJ  y/  =  O  ;  - (7) 

wAt/K,  va=  4-  -p  d  j^r  +  ^  /^z*" 
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It  is  convenient  at  this  point  to  give  the  solution 
where  a  periodic  vertical  force  with  magnitude  ZJ0(kr)  per  unit 
area  acts  along  the  interface  between  the  two  media.  The  functions 


Cp  and  must  satisfy  the  boundary  conditions: 

(>«.),  =  O  at  z  -  O - (*) 

(Prz)a  -  (Pzz),  =  ZJjKr)  aAr.-H - (9) 

(fz  r)a  -  O  O-t  Z  =  l-l  - (lo) 

Ur,  -  at  -2L-  H - (n) 


Typical  solutions  of  (5)»  (6)  and  (7)  are  of  the  form: 


Ja  (*r)  - - (/A) 

<Pa  r  Be~yZ  Jo  (Kr)  - - - - (13) 

VCe'rt  Ja(Kr), - - 


where  the  separation  constants  £  ,  T|  ,  J  ,  obtained  by  substituting 
(!2)-(l4)  in  (5)-(7)»  are 


8, 


£  ,  ly  ,  Y  'being  positive  real  or  positive  imaginary.  The 

time  factor  is  temporarily  omitted.  The  boundary  condition 

expressed  by  (8)  has  been  satisfied  by  the  choice  of  ^  .  If 

equations  (l2)-(l4)  are  substituted  in  the  boundary  conditions 
given  by  (9)-(ll)  using  (3)-(*0»  three  equations  for  the  constants 
A,  B,  C  are  obtained.  Solving  for  A,  B,  and  C  we  have  after  some 


reduction: 


/?  =  -= 


Z  CO 


V 


Mi.  4a  £  F  ( i<)  cotl  (g  h) 


B-  ^ 


Mi  Tf) 


(lb) 


C--  2 

Mi 


2.V  g 


v;here  F(K)  is  given  by: 


gM  =  A«(in)  -  -  (*  )  ] . (/7^ 


For  a  source  of  compressional  waves  beneath  the  inter¬ 
face  at  the  point  r  =  0,  z=H+d,  the  simplest  solution  which 
gives  a  compressional  wave  travelling  radially  outward  in  an  un¬ 


limited  solid  is 


i 
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where  R  =  ^  ^  denotes  the  radial  distance  from 

the  point  source.  If  an  image  is  placed  at  r  =  0,  z  =  H  -  ‘d,  we 
obtain  a  solution  for  which  the  tangential  stress  at  z  =  H  vanishes, 
namely 


-c 


Ml  o 


= 


(Jj  J 

.-'S'/ 


R 


•+ 


R‘ 


-  -  r/?; 


o 


where  R  =  ^ X*+y*+(z.  —  H  +c/)a  .  From  the  integral  form  of  the 
Bessel  function  Lamb^  showed  that  the  above  expression  for 


is  equivalent  to 


%  =  /  -  - 


V 

00 


JjKr)KJl<  +  r  e~V^z  ,  “ _  ^  (xr)KcJ, K 


=  zf(- e-’^r) kJk  : 


- (io) 


From  (3)  and  (4)  we  can  write  for  z  =  H 

fir  =0 


00 


Pzz  ^  c  ^ J0(Kr)K^K\ 


■(u) 


Referring  to  (9)  we  note  that  the  condition  for  conti¬ 
nuity  of  normal  stress  at  the  interface  will  be  satisfied  by 
superposing  on  (20)  the  solution  obtained  from  (l6)  and  (l2)-(l4) 
by  putting 


12.  Horace  Lamb,  loc.  cit. 


r 
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and  Integrating  with  respect  to  k  between  the  limits  0  and  oo  , 

It  can  he  verified  from  (2),  (12)-(l4)  and  (l6)  that  the  condition 
for  continuity  of  w  across  the  interface  will  he  satisfied.  Thus 
we  may  write  for  the  formal,  steady  state  solution  for  a  periodic 
compressional  point  source: 


. .(a «  -  c^J(Kr)  kJk 


00 

J  y 

o 


L 


-"-(43) 


(M) 


(iS) 
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Prom  (2),  (24)  and  (25)  we  may  write  for  the  displacements  q.g 
at  the  surface  of  the  solid  bottom  (z  =  H): 

00 


jF(K) 


Xa*M)  _  arje^j  ficr) k  d 

£4  ?  J 


< 


WH 


Integrals  of  this  general  type  have  been  evaluated  by 

13  14  15  l6 

Lamb  ,  Nakano  ,  Sezawa  and  Pekeris  .  The  procedure  is  to 

transform  the  path  of  integration  to  the  complex  k-plane.  The 

solutions  can  then  be  expressed  as  the  sum  of  the  residues  of  the 

integrands  and  two  integrals  along  branch  lines  corresponding  to 

the  branch  points  k  =  t*)/oc2  and  k  =  •  The  residues,  which 

diminish  as  l/^j"?  give  the  normal  mode  solutions,  whereas  the 

p 

branch  line  integrals  diminish  as  l/r  ,  and  become  negligible  for 
large  r. 

Following  this  procedure  and  neglecting  for  the  while 
the  contributions  of  the  branch  line  integrals,  we  find  for  large 
values  of  r: 


13.  Horace  Lamb,  loc.  cit . 

14.  H.  Nakano,  On  Rayleigh  Wave,  Jap.  Joum.  Astron.  and  Geophys.  , 

Vol.  II,  No.  5,  PP  233-326,  1925 

15.  KatsutB&a  Sezawa,  Love-Waves  Generated  from  a  Source  of  a 

Certain  Depth,  Bull.  Earth.  Res.  Inst.  Tokyo,  Vol.  13» 

Pt.  1,  pp  1-17 »  1935 

16.  C.  L.  Pekeris,  loc.  cit . 
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4.  n  \l7jr* r 


a 


lf(±>  (  4  4/  f, 


ft  at  £a*U„u)  ny 

a  t-  **  frt 


s  ^ 


-\d  i(at-K.„r-vfy) 

e  e  ....  (**) 


ur  - 

H  - 


an 


ur 

/t 


a_ 

7 rxnr 


^f^n) 


(5<n- 


^/4)eVnd  e(ujT~  Knr'm  —  -M 


The  subscript  n  indicates  that  the  quantity  is  to  he 
evaluated  at  k  =  kn,  where  kn  are  the  roots  of  the  period  equation: 

f(k)  -  o  - (3o) 

l  (wt-  Krr  - 11  fo) 

From  the  factor  £  it  can  be  seen 

that  the  phase  velocity  cR  can  be  expressed  by 

C  n  —  uu/ At n  _ _ _ _ _ '  -  (ilj 

and  (28),  (29)  can  be  written  in  the  form 

<?<K)  e’"‘'e,Y"t-',"r+v’)  . -W 

h  1/  /<  va 

"S'f  wWe-vVe'r‘rf''"r‘¥’)  . . w 

where 

^  =  (T^)  %  Cl  l4 

- H 


13. 


R.M  /3* 


(s- £/£) 


and 


rW= 


<?<£ 

*/£ 


r.  0-c'/«t)i  _  (*» H ) (/- Cn/< )vc« l 

_  4 


(Si) 


Equation  (30)  takes  the  dimensionless  form: 


4, [/-<■*/<  p-  °"A1  -  ( ■ a ' 


This  last  equation  defines  an  implicit  relationship 
between  the  period  "77  -  2/h and  the  phase  velocity 
cn,  with  the  elastic  constants  of  the  system  as  parameters.  It 
is  evident  that  the  wave  system  given  by  (28)  and  (29)  will  be 
attenuated  for  kn  complex,  the  degree  of  attenuation  depending 
on  the  magnitude  of  the  imaginary  component.  The  energy  loss 
associated  with  damped  propagation  is  due  to  the  refraction  of 
waves  from  the  liquid  layer  into  the  solid  bottom.  For  the  case 
>•  V)  which  is  considered  in  this  paper,  un¬ 

damped  propagation  will  occur  for  the  cases: 


X.  >  4z.  ^  C  ^  V, 

X~.  °<a  >  /SA  >  V,  >  C 
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Only  case  I  will  be  considered  in  this  paper. 

numerical  calculation  of  the  period  from  (37)  shows  it 


to  he  a  multiple  valued  function  of  phase  velocity,  each  value 
belonging  to  a  distinct  mode  of  propagation.  The  summations  in¬ 
dicated  in  (2g)  and  (29)  are  to  be  carried  over  all  the  modes 
n  =  1  to  n  =  cd  .  The  amplitude  factors  Q(kn)  and  W(kn)  give  the 
relative  excitation  of  each  mode  as  a  function  of  the  period.  The 
factor  €  ”  shows  the  influence  of  the  depth  of  the  source  on 
the  amplitudes  of  the  different  modes. 


It  can  be  shown^7  that  the  period  equation  (37)  expresses 


the  condition  of  constructive  interference  between  plane  waves 
undergoing  multiple  reflection  in  the  liquid  layer  at  angles  of  in¬ 


cidence  beyond  the  critical  angle  of  reflection  Oc  ,  where 

*  With  this  simplified  point  of  view  each 


mode  represents  a  different  order  of  interference  and  the  disturb¬ 
ance  at  a  distant  point  is  obtained  by  the  superposition  of  waves 
arriving  at  the  point  along  the  oblique  ray  paths  (defined  by  the 
angle  of  incidence  0  )  for -which  constructive  interference  occurs. 
For  the  case  of  plane  waves  the  wave  number  k  in  the  previous  dis¬ 
cussion  can  be  expressed  by  k  =  ,  where  1  is  the  wave¬ 

length  measured  along  an  oblique  ray  path.  The  plane  wave  theory 
however  fails  to  provide  an  adequate  explanation  of  how  waves  from 
a  source  in  the  bottom  can  enter  the  liquid  at  angles  of  incidence 
beyond  &c  • 


17.  Frank  Press  and  Maurice  Ewing,  loc.  cit. 
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Generalization  for  an  Arbitrary  Pulse 

In  an  dispersive  medium  in  which  an  arbitrary  initial 
disturbance  occurs,  the  energy  associated  with  each  period  is 
known  to  propagate  with  the  group  velocity  given  by  the  well  known 
formula 

c.  - - (}>) 

The  period  equation  was  used  to  obtain  the  phase  velocity 
as  a  function  of  kh.  The  group  velocity  was  subsequently  deter¬ 
mined  by  numerical  differentiation  according  to  (38).  The  results 
of  these  computations  for  c  ^  v^_  (case  I  above)  are  shown  in  Figures 
(2  )  and  (3  )  where  c/vq  and  U/v-^  are  plotted  as  functions  of  the 
dimensionless  parameter  S'-  tt(Jl  *  HfvtT  =  for  the 

cases  4//>  =4.$,**. /570*,  a*.  O**  V*  -  °<a.  -  J~3 . , 

@3LS  3  V)  which  represent  respectively  the  approximate  conditions 

for  a  granitic  and  basaltic  bottom.  Since  the  functional  dependence 
of  kh  on  c/v-^  is  not  very  sensitive  to  small  changes  in  the  value 
of  A/j?  we  have  taken  a  single  value  of  this  constant  for  both 

1® 

cases.  For  both  cases  two  modes  are  plotted.  In  a  previous  paper  0 
the  first  and  second  mode  was  considered  to  represent  ti^o  distinct 
branches  of  the  period  equation,  the  second  branch  having  higher 
modes.  This  distinction  is  artificial  and  unnecessary.  The  phase 
and  group  velocity  of  the  first  mode  (n  =  l)  approach  the  velocity  (Cr) 
of  Rayleigh  waves  in  the  bottom  layer  as  )>  — ^  O  ,  or  as  the  wave¬ 
length  becomes  very  long  in  comparison  to  the  thickness  of  the  first 
layer.  As  c/v^  and  U/v^— >  1,  Y  approaches  a  limiting  value.  For 

18.  Frank  Press  and  Maurice  Swing,  loc.  cit . 
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c/v^<  1  (case  II  above)  the  propagation  reduces  to  that  of  a  boundary 
wave  (Stoneley  wave)  along  the  bottom  interface.  In  the  second  mode 
(n  -  2)  the  phase  and  group  velocity  equals  the  velocity  of  shear  waves 
in  the  bottom  (c/v^  =  TJ/v-^  =  2  or  c/v-^  =  TJ/v-^  =  3  as  the  case  may  be)  for 
a  limiting  value  of  X  corresponding  to  a  cut-off  period  above  which  k  is 
complex  and  damping  occurs.  For  c  and  U  close  to  v^  the  periods  become 
infinitely  small.  Higher  modes  of  propagation  (n  =3»  ^  •••)  exist,  each 
having  the  same  cut-off  velocity  c  =  U  =  ,  but  increasingly  shorter  cut¬ 

off  periods.  In  general  the  periods  of  the  higher  modes  corresponding  to 
a  given  phase  or  group  velocity  become  progressively  shorter.  The  groupo 
velocity  of  the  first  and  second  modes  are  further  characterized  by  the 
existence  of  stationary  values  of  group  velocity,  the  importance  of  which 
will  be  discussed  in  a  later  paragraph. 

Having  obtained  the  solution  for  the  steady  state  propagation  of 
simple  harmonic  compressional  waves  originating  in  a  point  source  it  re¬ 
mains  yet  to  obtain  the  solution  for  an  arbitrary  initial  disturbance. 

If  the  time  variation  of  the  initial  disturbance  (19).  at  the  source  is 
f(t)  represented  by  its  Fourier  transform 


co 


-p(t)“  I  eCU)*  q_(u>)  c/uj 


00 


the  compound  displacements  q’  and  w’  for  an  arbitrary  pulse  can  be 

H  H 

written  as 


17. 


where  is  a  function  of  U>  through  (30)  and  g{td)  is  given  by 


-00 


19 

Pekeris  was  primarily  interested  in  a  source  consisting  of  an 

*  ~  cr  £ 

exploding  charge  and  he  appropriately  chose  -fvt)  =  £  for 

t  >  0  and  f(t)  =  0  for  t  <  0.  Our  primary  interest  is  the  impulse 
developed  in  an  earthquake  and  for  lack  of  more  detailed  informa¬ 
tion  we  choose  f(t)  to  vanish  for  all  but  infinitesimal  values  of 
t  in  such  a  manner  that  f (ft)  Jt  =  A  .  With  this 


choice  equations  (40)  and  (4l)  are  given  by: 


Q(k») 


CO - 


To  evaluate  these  integrals  we  follow  Pekeris  in  usin^ 
Kelvin's  approximate  method  of  stationary  phase.  The  method  assumes 
that  the  exponential  term  oscillates  rapidly  over  the  range  of  inte¬ 
gration,  whereas  the  remainder  of  the  intergrand  is  a  slowly 


19.  C.  L.  Pekeris,  loc .  cit . 


-  fra) 


-  -  fr3) 


-  -  frv) 


varying  function  of  v>.  The  principle  contribution  to  the  integral 
therefore  occurs  at  the  points  for  which  fa)ldu)  =  o  , 

where  "ffo)  .  An  expansion  of  f(w>)  about 

these  points  of  stationary  phase  is  made  and  formulae  for  the  approxi¬ 
mate  evaluation  of  the  integral  are  obtained.  Jeffreys^,  Lamb^, 

22 

and  Pekeris  derive  and  discuss  the  appropriate  formulae.  In 
evaluating  the  integrals  (I4.3)  and  (44)  we  shall  make  use  of  the 
equations  and  computations  as  given  by  Pekeris.  The  first  of  these 


can  be  written  as  follows: 

<50 

-  K(u>)r 


f  -  JrT3«i  VM  e 

-a,  I  dw*  1 


i(co0i-K(u>,)r  ifl/q) 


where  the  upper  or  lower  sign  is  to  be  taken  in  the  exponential 
according  as  is  positive  or  negative  and  k,  U,  0/ (jfc/ Y ,  C0o 

are  evaluated  from  the  group  and  phase  velocity  curve  for  the  values 

of  r  and  t  which  satisfy  r/ t  =  TJ.  In  deriving  the  second  equation 

_  y.  .  v  00  ^  I  I  -  du) 

of  (45)  use  was  made  of  the  relations  0  -  jjfv,  )  ^  ~  d~K  * 

U)  -  C  K  .  The  above  formula  is  valid  only  when  the  following 

condition  is  met: 


20.  Harold  Jeffreys,  loc.  cit . 

21.  Horace  Lamb,  Hydrodynamics,  Cambridge  University  Press,  6th 
Edition,  pp  395-398,  1932 

C.  L.  Pekeris,  loc.  cit. 


22. 


19. 
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c^Z 

dd-a 
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where 


z  - 


l/,a  i^K 

H  W1- 


JU 


XTtU3-  J  * 


7, 


This  condition  is  violated  for  small  r  and  for  points  near  a  sta¬ 
tionary  value  of  group  velocity.  In  the  latter  case  we  use  an  approx¬ 


imate  formula  devised  by  Airy: 


i//  \  c  (c0ot 

£M  Yfa)  e 


n°  ..  .  ,  ,  ,  **  £(v)  Y(&.)  e 
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E(v)  =  V,/}  [_J.‘/3(v)  +  J,h  M  J 
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^  U  vnj'o 

1  >  Tp 
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-t  >  -77 
t  <  H 

UvnA.% 
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-  (47) 


-  -  -  M 


--  (71) 
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Only  the  real  part  of  the  integral  in  (48)  has  been  considered. 

In  these  equations  the  subscript  o  indicates  that  the  quantity  to 

which  it  is  attached  is  to  be  evaluated  at  a  stationary  value  (min  or  max)  of 

group  velocity.  The  use  of  (4$)  is  justified  only  when 


where 

a.)  = 


^  c  u»,„ 


t 


r 

>  v. 


max. 


We  now  apply  these  results  to  equations  (43)  and  (44). 
In  order  to  compute  the  displacements  q^  and  w^,  one  first  finds 
the  contribution  of  the  hth  mode  at  a  given  time  t  and  range  r 
by  determining  ft),  ,  c/U-p  kH,  M r,/jr  from  the  phase  and  group 
velocity  curves  appropriate  to  the  medium  such  that  r/ t  =  U. 


21. 


Si 


*  and  Wg  are  then  obtained  by  summing  over  all  modes: 
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provided  r  is  large  and  t  is  sufficiently  removed  from  the  value 
t0  =  r/U0 ,  where  U0  denotes  a  stationary  value  of  group  velocity. 
We  have  added  to  the  expressions  (U3)  and  (44)  their  complex 
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conjugates  in  deriving  (56)  -  (59)»  since  Q  and  W  are  even 

functions  of  O),  f(u>)  is  an  odd  function  of  w*  and  there  are  two 

stationary  points  iu)  where  3Y^)/ du)*  o  # 

Near  the  time  t0  corresponding  to  a  stationary  value  of 

group  velocity  (Airy  phase)  one  first  determines  ,  fit),  k,  c,  Z, 

dz/d*  at  this  point  from  the  phase  and  group  velocity  curves 

1  1 

appropriate  to  the  medium.  As  before,  and  w0  are  then  computed 
hy  summing  over  all  modes: 


i  = 


uv,  A 


3Vs  H1  r 


<  Q_ (jA  M (MLfujt-Knr-hll/t/) 
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provided  r  is  lhrge  and  t  sufficiently  close  to  tQ  for  condition 
(53)  to  hold.  The  factors  E(v)  from  which  the  envelope  of  the 
Airy  phase  can  be  determined  and  G(v)/E(v)  from  which  the  range  of 

validity  of  (60)»  and  (6l)  can  be  found  was  plotted  by  Pekeris  and 

♦ 

is  reproduced  here  as  Figure  4. 
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Equations  (56)  through  (6l)  have  been  used  to  compute 
the  amplitudes  of  the  component  displacements  at  the  bottom  inter¬ 
face  for  an  impulsive  point  source  of  compress ional  waves  located 
at  a  depth  beneath  the  interface  equal  to  the  depths  of  the  water 
(d  =  H)  and  at  distances  of  1000  km.  and  10,000  km.  respectively. 

The  results  for  the  first  and  second  modes  are  presented  in 
Figures  5  to  8  for  the  case 

%  =  °(a.=f3/6a  83 ' 3V<  >  H  =  *  ;  r=  ;  /°.o OO 

n 

In  these  figures  the  amplitudes  are  plotted  as  functions  of  the 

—  "t  # 

dimensionless  parameter  — ^  where  t  denotes  the  time  after 

the  initial  impulse  and  is  the  travel  time  for  shear  waves 
in  the  bottom.  The  travel  time  of  waves  of  any  given  period  in 
the  first  two  modes  can  be  computed  from  Figure  9* 

In  order  to  show  the  effect  of  the  depth  of  the  source 
on  the  displacements  q^  and  w^  the  amplitudes  were  computed  for 
a  source  located  in  the  bottom  layer  close  to  the  interface  (d  =  0). 
The  resultant  amplitudes  in  the  first  two  modes  for  the  two  cases 
d  =  0  and  d  =  H  are  plotted  against  V--  w/i  in  Figures  10  to  13. 

Discussion  of  the  Normal  Mode  Solutions 

In  an  earlier  section  the  formal  steady  state  solutions 
(26)  and  (27)  were  expressed  as  the  sum  of  the  residues  of  the 
integrands  and  two  integrals  along  branch  lines  corresponding  to 
the  branch  points  c  =  oL^  and  c  =  .  The  residues  lead  to 


the  normal  mode  solutions  which  predominate  at  large  distances 


hecause  they  diminish  only  as  l/\[r  •  Now  the  normal  mode  solu¬ 
tions  vanish  at  c  =  U  »  0^  and  do  not  exist  for  c  >  (3%  .  In 
this  region  the  branch  line  integrals  contribute  waves  travel¬ 
ling  with  the  speed  of  compressional  and  shear  waves  in  the  bot¬ 
tom,  diminishing  rapidly  in  amplitude  with  distance  as  '/r*. 

At  a  time  "t  >  after  the  initial  impulse 

(i.e.,  immediately  after  the  arrival  of  shear  waves)  the  normal 
mode  solutions  predominate  providing  r  is  large.  Thus  the 
ground  motion  in  the  first  two  modes  due  to  a  distant  impulsive 
point  source  of  compressional  waves  having  a  moderately  broad 
spectrum  and  located  within  the  bottom  can  be  obtained  with  the 
aid  of  Figures  3  and  5-9* 

In  the  first  mode  the  group  velocity  approaches  the 
velocity  of  Rayleigh  waves  (  °r  =  •919^/S«  )  as  if” >0  •  The 
first  arrival  consists  of  long  period  Rayleigh  waves  whose  am¬ 
plitudes  and  period  vary  with  time  according  to  the  graphs  of 
Figures  5,  6  and  9*  At  the  time  t  =  £  a  short  period  wave 
travelling  with  the  speed  of  sound  in  water  arrives*.  According 
to  Figures  10,  11  amplitudes  associated  with  this  wave  are  par¬ 
ticularly  small  unless  the  source  is  located  close  to  the  inter¬ 
face.  This  follows  from  the  factor  in  equations  (56)-(6l) 

which  gives  the  dependence  of  the  normal  mode  amplitudes  on  the 


*  The  group  velocity  curves  show  an  upper  frequency  limit  for 
these  waves.  For  higher  frequencies  the  solutions  reduce  to 
boundary  waves  propagated  along  the  interface  in  a  manner 
analagous  to  Stoneley  waves. 


25. 


depth  of  the  source  "beneath  the  interface.  Since 


the  amplitudes  of  the  shorter  period  waves  (  large)  are  cut 
down  if  the  source  has  any  depth  x^hatever.  For  source  depths 
greater  than  several  times  the  depth  of  water  the  amplitudes 
of  all  "but  the  longest  period  waves  are  reduced  to  insignifi¬ 
cantly  small  values. 


the  short  period  and  long  period 


For 


"branches  of  the  group  velocity  curve  approach  each  other  and 
finally  merge  at  a  minimum  value.  The  Airy  phase  associated 
with  this  stationary  value  of  group  velocity  consists  of  waves 
of  large  amplitude  having  a  prescribed  period  and  velocity.  Since 
waves  propagated  at  a  stationary  value  of  group  velocity  decrease 
with  distance  as  /  f  rS^  (equations  60  and  6l)  and  waves  pro¬ 
pagated  at  other  values  of  group  velocity  diminish  as  l/r,  the 
relative  amplitude  of  the  Airy  phase  increases  with  distance. 

This  effect  is  shown  in  Figures  5  and  6  where  the  amplitudes  have 
been  computed  for  the  two  ranges  r  =•  1000  km.  and  r  =  10,000  km. 

The  ground  motion  in  the  second  mode  (Figures  J,  g,  12, 
13)  begins  with  waves  arriving  with  a  cut-off  period  at  a  speed 
(J  =  ^.  The  amplitudes  are  zero  at  the  onset  and  thereafter  in¬ 
crease.  At  a  time  corresponding  to  propagation  at  the  maximum 
value  of  group  velocity*  large  amplitude  waves  arrive  with  an 
abrupt  onset.  Thereafter  these  waves  separate  into  several 


*  Pekeris  attaches  the  name  Airy  phase  in  reverse  to  maximum 
values  of  group  velocity. 
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simultaneous  arrivals  of  different  periods  and  amplitudes.  At 
the  time  t  =  r/v^  a  short  period  wave  (  JT-*  oo  )  arrives  travel¬ 
ling  with  the  speed  of  sound  in  water.  As  was  the  case  in  the 
first  mode,  the  amplitude  of  this  phase  is  negligibly  small 
unless  the  source  is  located  at  or  very  near  to  the  interface. 

For  t  >  r/v^  the  two  arrivals  corresponding  to  two  branches  of 
the  group  velocity  curve  approach  each  other  and  merge  at  a 
minimum  value  of  group  velocity,  producing  an  additional  Airy 
phase  in  the  second  mode.  Since  the  second  mode  Airy  phase 
occurs  at  a  shorter  period,  its  amplitude  is  smaller  than  that 
of  the  first  mode  if  the  source  depth  is  of  the  order  of  the 
thickness  of  the  liquid  layer  or  greater. 

Table  I  gives  the  period  T  and  group  velocity  TJ/v^ 
associated  with  each  of  the  stationary  values  of  group  velocity 

of  the  first  two  modes  for  the  two  cases  a)  =  2  »  S 

=  l/Tfa l;  fa  -  2  V,  y  V,  -  JfClo  }  H  =  S /* , 

and  b)  Pijp, =  U  =  Sozx> ***£**  * 


TABLE  I 


Periods  and  G-roup  Velocities 

of  Airy  Phases 

y 

T 

u/v, 

Case  a  First  Mode  .33 

10.  sec. 

.SO 

Second  Mode  *49 

6.7 

1.7 

3.^ 

■75 

Case  b  First  Mode  .  2S 

12 

.64 

Second  Mode  .46 

7.2 

2.6 

.82 

4.0 

.56 
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The  minimum  value  of  group  velocity  in  the  second  mode  indidated 
by  the  dashed  lines  in  Figures  2  and  3  is  not  included  in  Table  I 
since  its  exact  position  has  not  been  calculated.  We  infer  that 
amplitudes  associated  with  this  point  are  not  significantly  large 
because  of  its  proximity  to  the  cut-off  period  and  its  extremely 
sharp  curvature. 

The  ratio  of  vertical  to  horizontal  displacement  for  the 
first  two  modes  is  plotted  as  a  function  of  in  Figure  14.  It 
is  seen  that  for  o  the  ratio  -y-  approaches  the  correct 

value  for  Rayleigh  waves  propagated  along  the  surface  of  a  semi¬ 
infinite  elastic  solid  for  which  )\-  JU  . 

From  equations  (2),  (3)»  and  (23)  it  can  be  seen  that 
the  factor  sin  (£^2-  )  gives  the  vertical  variation  of  pressure 
and  horizontal  displacement  and  cos  (^2: )  gives  the  vertical  vari¬ 
ation  of  vertical  displacement  in  the  liquid  layer.  Nodes  of 
pressure  correspond  to  antinodes  of  vertical  displacement  and  vice 
versa.  In  Figure  15»  the  vertical  pressure  variation  in  the  liquid 
layer  is  presented  as  a  function  of  X  for  the  first  two  modes. 

The  ground  motion  at  a  point  is  obtained  by  the  super¬ 
position  of  the  contributions  of  all  modes.  Our  discussion  thus 
far  has  been  limited  to  the  first  two  modes.  The  periods  of  the 
higher  modes  corresponding  to  a  given  group  velocity  become  pro¬ 
gressively  shorter.  For  source  depths  of  the  order  of  a  fraction 
of  the  depth  of  water  or  more,  the  higher  modes  will  contribute 


little  to  the  ground  motion  because  of  the  factor  £ 


In 


addition,  the  contributions  of  the  short  period  higher  modes  will 
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"be  cut  down  because  of  their  greater  susceptibility  to  scattering 
by  an  irregular  ocean  bottom.  The  poor  short  period  response 
characteristics  of  intermediate  and  long  period  seismographs  will 
further  discrimate  against  the  higher  modes. 

The  theory  and  computations  described  represent  an 
earthquake  by  an  impulsive  point  source  of  compressional  waves. 
Actually  an  earthquake  is  neither  of  these.  A  true  source  is 
usually  extended  and  rich  in  distortional  waves. 

Despite  the  highly  idealized  assumptions  concerned  the 
nature  of  an  earthquake  source,  it  is  believed  that  the  normal 
mode  theory  given  above  is  sufficiently  adequate  to  explain  the 
salient  features  on  the  seismograms  of  shallow  focus  submarine 
earthquakes. 

Thus  far  only  the  Airy  phase  of  the  first  mode  and  the 

short  period  water  wave  phase  (MT'*  phase)  have  been  observed.  In 

the  next  section  the  data  on  the  Airy  phase  is  presented.  The 

2"z 

"TM  phase  has  been  described  in  an  earlier  report.  J 


23.  The  MTM  Phase  of  Shallow  Focus  Submarine  Earthquakes,  Technical 
Report  No.  1  by  Ivan  Tolstoy,  Maurice  Ewing  and  Frank  Press. 
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The  Airy  Phase  Observed  at  Bermuda  and  Euancayo 
1.  Bermuda. 

The  Bermuda  seismological  station  is  situated  at  32*23  *N, 
64*4l*W.  It  is  equipped  with  a  Milne-Shaw,  two  components, 
oriented  NE-SW  and  NW-SE.  The  free  period  is  given  as  12  sec  and 
the  magnification  as  250  by  the  USCGS  seismological  bulletin. 

A  series  of  records  was  examined  in  the  hope  of  detecting 
late  arrivals  of  9  -  12  sec  period.  Such  arrivals  were  detected 
on  all  readable  seismograms  of  the  Dominican  Republic  shocks,  but 
usually  not  for  shocks  occuring  elsewhere  in  the  Atlantic.  Thus 
a  series  of  Mid-Atlantic  Ridge  quakes  were  studied,  including  the 
large  shocks  of  April  11,  1946  (leS,  l4  l/2*W,  Pasadena  magnitude 
7.2)  and  November  25,  19^1  (37  l/2*N,  IS  l/2*W,  magn.  8  l/4),  and 
with  one  exception,  all  gave  negative  results.  We  will  dwell  upon 
this  point  later. 

Table  II  gives  the  results  obtained  for  a  series  ot 
Dominican  Republic  shocks.  H  is  the  origin  time,  Ma  gives  the 
arrival  time  of  the  maximum  amplitude  wave  of  the  group  inter¬ 
preted  as  being  the  Airy  phase.  It  will  be  noted  that  in  some 
cases  these  times  appear  to  differ  somewhat  for  the  two  components. 
The  column  marked  t.t  gives  the  total  travel  time  of  the  Airy 
phase  in  minutes  and  seconds.  Aa  gives  the  maximum  double  trace 
amplitude  of  the  Airy  phase,  Ap  and  As  give  those  of  the  P  and  S 
groups  (in  cm).  T  gives  the  period  of  the  group  on  the  seismogram 
in  seconds.  This,  as  well  as  the  arrival  times  are  only  approxi¬ 
mate,  since  the  paper  velocity  of  these  seismograms  is  only  Smm/min. 
However  they  are  sufficient  for  the  present  purposes. 
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This  list  includes  only  shocks  of  magnitude  7  or  less. 

The  large  shocks,  such  as  those  of  Aug  8,  1946,  (magn.  7*6)  and 
Aug  4,  1946,  (magn.  8),  July  29,  1943»  (magn.  73)  and.  others  are 
not  included.  The  reason  for  this  is  the  impossibility  of  distin¬ 
guishing  the  phase  we  are  studying,  due  to  the  very  large  trace 
amplitudes  and  consequent  confusion  and  fading  of  the  traces. 

Arrival  times  and  velocities. 

Let  us  see  now  what  these  arrival  times  correspond  to. 

The  times  given  in  Table  II  were  read  for  the  maximum  amplitude  of 
the  group  of  waves  interpreted  as  Airy  phases.  Such  a  procedure  is 
justified  by  the  theory  of  wave  propagation  in  dispersive  media:  for 
a  stationary  value  of  group  velocity  there  is  a  gradual  build  up  of 
the  amplitudes  until  the  arrival  of  the  corresponding  wave  length, 
at  which  point  the  amplitude  is  maximum. 

For  purposes  of  comparison  of  the  arrival  times  of  this 
phase  with  those  of  the  classical  earthquake  phases,  it  would  be 
necessary  to  pick  the  time  at  the  beginning  of  the  group.  For  this 
one  should  subtract  roughly  2  min  from  the  preceeding  arrival  times. 

The  first  remark  that  we  shall  make  is  that,  even  with  the 
latter  provision,  the  travel  times  are  much  too  large  to  be  any  of 
the  known  tabulated  earthquake  phases^^  such  as  either  multiply 
reflected  body  waves,  Love  and  Rayleigh  waves.  They  are  too  large 
for  ScS  even  if  this  could  be  observed  at  such  a  short  distance  from 
the  epicenter. 


(24) 


s.and  K.  E,  Bullen,  Seismological  Tables,  British 
for  the  Advancement  of  Science,  1940 
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Also  these  phases  are  not  aftershocks,  as  could  easily  be 

verified  by  examining  the  records  obtained  at  sei sinological  stations 

having  instruments  with  different  characteristics  (such  as  the  Weston 

(25)  (26) 

and  Fordham  Benioffs).  As  the  reader  may  verify, /the  published 
lists  of  aftershocks  do  not  indicate  any  shocks  whose  arrival  times 
would  fit  those  of  the  Airy  Phase, 

We  purposely  did  not  include  the  aftershock  of  Aug  7»  19^6* 
H  =  19:21:32  in  our  list,  since  it  was  followed  by  another  shock  15 
min  35  sec  later  which  did  interfere  with  the  Airy  phase  of  the 
former. 

We  see  therefore  that  this  arrival  is  clearly  a  new  phase, 
which  cannot  be  explained  by  any  of  the  standard  types  of  elastic 
waves  transmitted  through  the  earth!s  outer  or  inner  layers. 

Figure  l6  a,b  reproduces  seismograms  for  both  components 
of  the  shocks  of  Aug  8,  19*+6,  (H  =  17:24:07)  and  Aug  9»  1946,  (H  = 
08:25*40).  The  phase  in  question  is  seen  to  stand  out  very  conspi¬ 
cuously  from  the  late  coda.  On  all  records  examined  it  stands  out 
with  the  same  perfect  dlarity,  and  even  for  shocks  much  smaller  than 
any  of  those  enumerated  it  is  possible  to  detect  it,  though  it  is 
not  easy  to  get  good  times. 

To  obtain  the  velocity  of  this  phase  we  should  have  a 
fairly  accurate  idea  of  the  distance  from  the  epicenter.  Unfor¬ 
tunately,  it  is  not  possible  to  obtain  the  accurate  readings  of  the 
P-S  interval  necessary  for  this  purpose  from  this  type  of  seismogram. 
We  will  have  to  be  content  with  an  approximate  evaluation. 


(25)  N^E.S.A,  Bulletin,  kug.  1946 

(26)  Bull.  Seism.  Lab.,  Cal.  Inst,  of  Tech.,  July-Sept,  19*46 
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The  largest  of  the  shocks  listed  in  Table  II  is  that  of 
Oct  4,  19^6,  H  =  l4:45:32.  Its  epicenter  is  given  by  Gutenberg  as 
18*45  *N  6g*30lW»  This  means  A  =  l4°03l.  Taking  the  travel  time 
as  2^min  31sec,  we  get  a  velocity  of  3480  feet/sec  which  is  0.7  of 
the  speed  of  sound  in  water. 

Considering  the  epicenter  to  have  been  determined  to  within 
-  1*,  which  is  a  reasonable  estimate,  the  above  velocity  is  accurate 
to  within  £  250  feet/sec,  i.e.  ,  it  is  somewhere  between  0.75  and 
O.65  of  the  velocity  of  sound  in  water* 

The  period  of  this  phase  being  on  the  record  of  the  order 
of  9  sec,  it  is  seen  that  it  is  in  approximate  agreement  with  the  theore¬ 
tically  predicted  Airy  phase.  It  might  be  necessary  to  assume  a 
slightly  higher  bottom  velocity  than  that  of  granite  between  Bermuda 
and  the  West  Indies,  in  order  to  explain  its  slightly  lower  velocity. 

As  for  the  differences  in  travel  times  of  the  Airy  phases 
of  Table  II,  it  should  be  remembered  that  they  are  sensitive  to 
differences  in  distance,  since  the  velocity  of  propagation  is  com¬ 
paratively  low.  Thus  a  difference  in  distance  from  the  station  to 
the  epicenters  of  two  shocks  of  1*  would  cause  a  difference  in  travel 
times  of  the  order  of  1  min  4o  sec.  This  is  of  the  order  of  magni¬ 
tude  of  the  differences  shown  in  Table  II. 

Bermuda  is  obviously  ideally  located  for  the  recording  of 
Airy  phases  of  the  Dominican  Republic  shocks,  since  the  entire  paths 
are  oceanic.  However  a  more  favorable  type  of  instrumentation  would 
be  highly  desirable  in  order  to  make  a  thorough  study  of  the  Airy 
phase.  In  particular  a  vertical  component  seismometer,  with  a  peak 
sensitivity  in  the  vicinity  of  llsec,  a  higher  magnification  and 
greater  paper  speed  would  be  useful  in  obtaining  more  accurate  data. 
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These  Airy  phases  do  not  show  up  on  the  Milne-Shaw  seismo¬ 
grams  of  the  Dominion  Observatory  in  Ottawa.  Probably  an  appreciable 
amount  of  the  energy  of  the  phase  is  lost  in  the  transfer  from  ocean 
to  continent. 

A  fair  sized  shock  occurred  on  the  western  flank  of  the 
Mid-Atlantic  Ridge  on  May  6,  1944,  at  22*  .4  N  44.g*  W  (USCGS) 

H  *  00:13*7.  This  appears  to  be  the  only  one  of  the  series  of  Mid- 
Atlantic  Ridge  quakes  in  recent  years  to  have  produced  a  recognizable 
Airy  phase  on  the  Bermuda  seismograms.  The  amplidude  of  the  group 
of  waves  interpreted  as  such  is  small;  nevertheless  this  group  is 
fairly  well  individualized  and  stands  out  well  from  the  rest  of  the 
late  coda.  The  period  is  9  sec,  and  the  arrival  time  00:46:27, 
which  means  a  travel  time  of  32  min  45  sec  for  A  =  20e  17%  i.e. , 
a  velocity  of  37^0  feet/sec,  0.75  of  the  velocity  of  sound  in  water. 
The  double  trace  amplitudes  in  cm.  refd^As  follows: 
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2 .  Huancayo 

The  Huancayo  Magnetic  Observatory  is  situated  at  12*03 rS, 
75<’20*f.  It  is  equipped  with  2  Wenner  seismographs  oriented  NS  and 
EW.  The  natural  period  is  given 'as  10  sec  and  the  magnification  as 
1550  by  the  USCGS  seismological  bulletin.  It  also  has  a  Benioff 
vertical  of  period  1  sec  and  magnification  unknown. 
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On  Aug  5.  1944 ,  both  components  of  the  Wenner  recorded 
what  appears  to  be  the  Airy  phase  of  a  shock  which  occured  at 
01:24:08  at  13  l/2cS  92  l/2*W  (USCGS),  south  west  of  the  Galapagos 
Is.  The  arrival  time  was  01:50:42,  which  gives  a  travel  time  of 
26  min  34  sec.  Since  A  *  17*  the  velocity  comes  out  to  be  3^90 
feet /sec  or  0.8  of  th'e  velocity  of  sound  in  water*  The  period  is 
8  sec.  Figure  17  shows  this  Airj;  phase.  It  is  seen  to  be  very 
conspicuous.  The  amplitudes  in  cm.  are: 
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It  did  not  record  on  the  Benioff  vertical,  the  period  of  which  is 
too  short.  This  latter  record  proves  conclusively  that  this  is 
not  a  separate  shock,  but  a  phase  of  the  quake  recorded  selectively 
only  on  the  Wenner. 

Discussion 

The  velocities  obtained  for  the  Airy  phase  at  Bermuda  for 
the  Oct  4,  1946,  Dominican  Republic  and  the  May  6,  1944,  Mid-Atlantic 
Ridge  quakes  and  at  Huancayo  for  the  Aug  5.  1944,  Galapagos  shock  are 
respectively  3480  *  250  feet/sec,  3760  and  3^90  feet/sec. 

The  second  and  third  shocks  were  smaller  than  the  first, 
so  that  in  spite  of  the  greater  paths  travelled  in  each,  the  estimated 
velocities  cannot  be  considered  as  being  any  more  accurate.  The  dif¬ 
ference  between  the  first  two  velocities,  3480  -  ?50  and  3760  is 

therefore  sufficiently  small  to  be  explainable  by  errors  in  the 
epicentral  determinations. 
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As  for  the  velocity  obtained  from  the  Galapagos  shock,  it 
should  be  remembered  that  Huancayo  is  situated  a  little  over  150 
miles  inland,  so  that  the  velocity  of  3^90  feet/sec  is  probably  a 
little  high.  It  appears  only  reasonable  to  assume  that  once  the  Airy 
phase  type  motion  reaches  land  it  starts  travelling  with  a  velocity 
similar  to  the  velocity  of  microseisms  on  land,  roughly  3  times  as 
fast  as  its  velocity  along  an  ocean  path.  Allowing  for  the  path 
acfoss  land,  the  velocity  of  the  Airy  phase  through  the  ocean  for 
this  shock  is  .7  of  the  velocity  of  sound  in  water.  There  are 
therefore  no  real  discrepancies  in  our  figures  for  the  velocities. 

As  for  the  amplitudes  of  the  various  Airy  phases,  it 
appears  that  when  the  epicenter  of  tne  shock  is  in  approximately 
the  same  azimuth  as  one  of  the  components  of  the  recording  instru¬ 
ment,  it  is  this  component  that  comes  out  largest.  This  is  the 
case  for  the  Mid-Atlantic  Ridge  quake  of  May  6,  1544,  (the  NW 
component  of  the  Airy  phase  is  the  largest),  and  for  the  Galapagos 
quake  (the  EW  component  is  predominant).  These  two  records  appear 
therefore  to  suggest  that  in  the  horizontal  plane  the  motion  is 
mostly  longitudinal.  Because  of  their  orientation  the  Bermuda 
records  do  not  yield  any  information  on  that  score  (Table  Ilshows 
that  the  amplitudes  are  about  the  same  on  both  components.) 

We  have  already  mentioned  the  fact  that  for  the  two  large 
shocks  of  Nov  25,  19^1,  (37  l/2*  »,  18  l/2e  W,  magn.  8  l/4)  and 
April  11,  1946,  (1®S,  l4  l/2°W,  magn.  7.2)  the  Bermuda  Milne-Shaws 
failed  to  record  any  Airy  phases.  A  bathymetric  chart  of  the 
Atlantic  Ridge.  We  believe  that  this  provides  an  adequate  explanation 
of  the  failure  of  the  Airy  phases  to  show  up  in  the  Bermuda  instruments. 
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The  Mid-Atlantic  Ridge  forms  a  continuous  barrier  across  the 
Atlantic,  100  or  more  miles  in  width,  rising  well  above  the  2000 
fathom  line  throughout  its  whole  length,  except  for  one  or  two  narrow 
trenches.  In  sections  it  rises  above  the  1000  fathom  line.  It  is 
to  be  expected  that  an  Airy  phase  reaching  this  barrier  would  lose 
much  of  its  energy.  What  would  be  deft  after  crossing  this  barrier 
would  not  have  enough  amplitude  to  reach  Bermuda. 

As  for  shocks  occur ing  on  the  western  flank  of  the  Ridge, 
all  of  those  which  occured  in  recent  years  would  have  been  too 
small  to  produce  an  Airy  phase  with  appreciable  amplitude. 

It  also  seems  probable  that  the  Airy  phase  does  not  pro¬ 
pagate  over  distances  of  the  order  of  several  thousand  miles,  or 
at  least,  if  it  does,  its  amplitude  is  not  large  enough  to  dis¬ 
tinguish  it  from  the  very  late  coda.  Thus  the  torsion  seismometers 
and  linear  strain  seismometers  of  the  west  coast  do  not  show  any 
Airy  phases  for  the  Alaskan  and  Japanese  shocks,  and  neither  do 
the  Honolulu  Milne  Shaws;  the  latter  it  is  true,  are  troubled  by 
a  very  high  background. 

Mechanism  of  production  of  an  Airy  phase. 

It  can  be  shown  that  an  impulse  either  in  the  liquid 
layer  or  in  the  bottom  will  produce  waves  that  will  propagate 
according  to  the  laws  of  normal  mode  theory. 


Thus  any  earthquake  whose  epicenter  is  at  sea  or  close 
to  the  shore  will  produce,  among  other  types  of  waves,  an  Airy 
phase  after  propagation  for  a  considerable  distance  through  the 
acoustic  system  consisting  of  ocean  and  bottom.  However  it  will 
be  only  under  certain  given  conditions  that  this  phase  will  have 
sufficient  energy  to  produce  a  noticeable  effect  on  a  seismometer 
situated  at  a  certain  distance. 

It  will  be  noted  that  all  the  shocks  for  which  we  hare 
found  Airy  phases  are  shallow  focus  ones.  The  focus  is  therefore 
situated  practically  on  the  liquid-solid  interface.  The  impulse  on 
the  liquid  layer  is  strongest  under  these  conditions.  It  is  probabl 
that  the  Airy  phase  produced  by  a  shock  whose  focal  depth  is  appre¬ 
ciably  over  5-10  km.  is  too  weak  to  be  recorded  at  any  distance. 

All  our  data  »*efer  to  records  obtained  by  seismometers 
situated  on  land.  It  seems  quite  probable  that  with  the  aid  of 
sensitive  hydrophones  in  the  ocean,  at  sufficient  depth  to  remove 
the  effect  of  surface  gravity  waves,  it  would  be  possible  to  record 
the  Airy  phases  of  a  great  many  more  shocks  and  at  considerably 
greater  distances. 

Practical  applications 

The  study  of  the  Airy  phase  of  a  sufficient  number  of  sea 
quakes  would  provide  interesting  information  concerning  the  nature 
of  the  ocean  bottom,  since  the  latter  has  an  influence  upon  both 
the  periods  and  the  velocities  of  the  Airy  phase. 
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Conversely,  once  the  velocities  of  the  Airy  phase  are 
known  for  certain  paths,  the  recording  of  these  hy  means  of  suitable 
instruments  could  be  useful  for  the  more  accurate  determination  of 
epicentral  distances  and  epicenter  locations. 


% 
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Figure  17.  Airy  Phase  recorded  at  Huancayo  for  the  earthquake  of  1944  Aug  5  17 
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